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We investigate the existence of spatially localised solutions, in the form of discrete breathers, in 
general damped and driven nonlinear lattice systems of coupled oscillators. Conditions for the ex- 
ponential decay of the difference between the maximal and minimal amplitudes of the oscillators are 
provided which proves that initial non-uniform spatial patterns representing breathers attain expo- 
nentially fast a spatially uniform state preventing the formation and/or preservation of any breather 
solution at all. Strikingly our results are generic in the sense that they hold for any, attractive in- 
^^ • teraction, coupling strength and on-site potential and general driving fields. Furthermore, our 

I ' rigorous quantitative results establish conditions under which discrete breathers in general damped 

("—^ , and driven nonlinear lattices can exist at all and open the way for further research on the emergent 

rvq ■ dynamical scenarios, in particular features of pattern formation, localisation and synchronisation, 

, ' in coupled cell networks. 

TpT. PACS numbers: 05.45.-a, 63.20.Pw, 45.05.-fx, 63.20.Ry 

p^ ■ Intrinsic localised modes (ILMs) or discrete breathers in nonlinear lattices have attracted significant interest recently, 

not least due to the important role they play in many physical realms where features of localisation in systems of 

I""! coupled oscillators are involved (for a review see [l| and references therein), 0- [J]. For conservative systems proofs 

^■^ , of existence and (exponential) stability of breathers, as spatially localised and time-periodically varying solutions, 

*^ ' were provided in Q and respectively. Analytical and numerical methods have been developed to continue breather 

Ph ', solutions in conservative and dissipative systems starting from the anti-integrable limit 7] , [8] . During recent years the 

^J^ ' existence of breathers has been verified in a number of experiments in various contexts including micro-mechanical 

^ C^ l cantilever arrays [9], arrays of coupled Josephson junctions (l(|. cout)led optical wave guides [IJL Bose-Einstein 

condensates in optical lattices [13], in coupled torsion pendula |14| . electrical transmission lines llSj, and granular 

r\i crystals [l6[. Regarding their creation mechanism in conservative systems, modulational instability (MI) provides 

^ . the route to the formation of breathers originating from an initially spatially homogeneous state imposed to (weak) 

^D ' perturbations. To be precise, the MI of band edge plane waves triggers an inherent instability leading to the formation 

C ■ [ of a spatially localised state [l3]- Departing from this often too ideal assumption of a conservative system requires, 

^~1^ ■ for more realistic models, the inclusion of dissipation. Accomplishing breather solutions in the presence of dissipation 

^^ ' requests some compensating energy injection mechanism. As far as their persistence is concerned it is expected that 

\l ■ breathers can be continued from a conservative system into a system augmented by weak dissipation and driving. 

^— ^ Compared to their Hamiltonian (conservative) counterparts breathers in dissipative and driven lattice systems do 

, , not occur in families of localised solutions as they are provided by discrete sets of attractors for appropriate initial 

conditions contained in the corresponding basin(s) of attraction |18|-[25|. 

The aim of this work is to establish quantitative conditions in parameter space for the existence respectively non- 
existence of discrete breathers in general damped and driven anharmonic lattice systems. To this end we show that 
there exist parameter ranges such that for any launched localised state the difference between the maximal amplitude 
d • and the minimal amplitude of the oscillators decays exponentially fast. Consequently, a spatially uniform state is 
attained. Most importantly, this rules out the persistence of any non-uniform pattern such as breathers. Moreover, 
our results also identify parameter ranges for which no inherent instability that is able to trigger the formation of a 
localised pattern exists. Crucially, our rigorous quantitative results establish prerequisites for the existence of discrete 
breathers in general damped and driven nonlinear lattices beyond the validity of the continuation process starting 
from the anti-continuum limit [8|,[18j-[2l|. Our results are generic as they hold not only for any coupling strength 
but also for any hard on-site potential, any type of attractive interaction, any degree of attractive interaction, general 
driving fields and any breather solution. 

We study the dynamics of general driven and damped nonlinear lattice chains given by the following system 

Qn = -U'{qn)-iqn + A{t) + B{t)q^ 

r 

- ^l^n[V'{qn+3-qn) + V'{qn-qn-])\- (1) 

i=i 

The variable qn{t), where n G Z, is the amplitude of the oscillator at site n evolving in an anharmonic on-site potential 
U{qn)- The prime ' stands for the derivative with respect to <?„ and an overdot 'represents the derivative with respect 



- ' — \ 



to time t. 

The on-site potential U is analytic and is assumed to have the following properties: 

C/(0) = t/'(0) = , C/"(0) > 0. (2) 

In what follows we differentiate between soft on-site potentials and hard on-site potentials. For the former (latter) 
the oscillation frequency of an oscillator moving in the on-site potential U{q) decreases (increases) with increasing 
oscillation amplitude. A soft potential possesses at least one inflection point. If a soft potential possesses a single 
inflection point, denoted by qi, we suppose without loss of generality (w.l.o.g.) that qi > 0. Then the following 
relations are valid 

U'{-oo <q<0) < 0, C/'(0 <q<q^)>0, (3) 

U"iq,) = 0, C/"(-^ <q<q,)>0. (4) 

If U{q) possesses two inflection points denoted by qi^^ < and qi^+ > it holds that 

f/'(g,,_ <q<0) < 0, [/'(O <q< g,,+) > (5) 

f/"(*,±) = 0, U"{q,,^ < q < g,,+) > 0. (6) 

We remark that U{q) can have more than two inflection points (an example is a periodic potential U{q) — — cos{q)). 
However, in the frame of the current study we are only interested in motion between the inflection points adjacent to 
the minimum of U{q) at g = 0. Hence, in the forthcoming we suppose that for soft on-site potentials the motion at 
each lattice site n is stays inbetween the inflection points, viz. q^,^ < qn(t) < qi^+, where U{q) is convex. 

Hard on-site potentials are, in addition to the assumptions in ^, characterised in their entire range of definition 

by 

U'{q < 0) < 0, U'{q > 0) > 0, U"{q) > 0. (7) 

In contrast to soft potentials, since the hard potentials are by assumption convex in their entire range of definition 
no boundedness condition as in relation to the motion in soft potentials is required. 

Each oscillator interacts to either side with its r neighbouring oscillators with (local) coupling strength k„ > (r 
determines the interaction radius) via forces derived from an attractive interaction potential V{u) which is analytic 
and furthermore, is assumed to have the following features: 

V{0) = V'{0) = 0, V"{0) > 0, V"{u / 0) > 0. (8) 

Thus V{u) is convex which is further characterised by V'{u > 0) > and V'{u < 0) < 0. It is through the site- 
dependent coupling strength that heterogeneity enters the model. The interaction potential can be harmonic but also 
anharmonic such as in /?— Fermi-Pasta-Ulam chains and Toda-type interactions. 

The parameter 7 > regulates the strength of the damping. A{t) and Bit) are smooth functions representing 
general external time-dependent fields with 

maxA(t) = Amax < 00, mm A{t) == Amtn > -00, 
tern teR 

maxi?(t) = Brnax < oo , mm B{t) = Bmin > -00. 

teR teR 

The A{t) and B{t) term in Eq. ([T]) is associated with direct and parametric driving respectively. 

For the sake of our study, that is to prove the absence of localised spatial patterns, we consider single breather 
solutions being peaked around one lattice site and the oscillators perform in-phase motion. Later we comment on 
the extension of our approach to multibreathers. That is, Eq. ([ij is supposed to possess a time-periodic and spatially 
exponentially localised solution, viz. a discrete breather, qn(t + Th) = (?„(i), with period T^ = 27r/wh where ojb denotes 
the breather frequency. Unless specified otherwise we refer in the following to such a solution as a breather. In 
general, breathers, being supported by periodic closed orbits in phase space, are associated with periodic bounded 
motion of the oscillators inside their on-site potentials U{q). Periodic solutions require time-periodic external fields 
A{t) = A{t + Ta) and B{t + Tg) with appropriate periods Ta and Tb- 

We introduce the following quantity related to the extremal values of the coordinates: 

qmax (t) = max g„ (i) , g„„;„ (i) = min (?„ (i) , 

n n 

and denote the difference between the corresponding velocities as 

Aq(t) = qmax{t) - qrnin{t) > . 



Non-uniform (uniform) states are characterised by non- vanishing (identicaUy vanishing) Aq{t). 
In the following we list the conditions satisfied by breather solutions described above: 

For breather solutions with period Tj, it holds that Aq{t) — Aq{t + Tb) and Aq{t) = Aq{t + Ti,). Furthermore, Aq{t) 
involves the same two oscillators all the time. In fact, since the breather pattern is spatially exponentially localised 
the two lattice sites involving qmaxit) and qmin(t) remain the same and only exchange their role after every change 
of sign of the periodically oscillating amplitudes. To be precise, the lattice sites n — max and n = min supporting 
the oscillators with qmax and q-min respectively during phases when (/n > swap when the coordinates qn become 
negative. To describe the behaviour of Aq{t) and Aq{t) we express a period duration Th = Td+ (71, — Td) as the sum 
of two stages of length Td and Tb — Td during which the coordinates qn{t) possess opposite sign. We first consider 
stages of length Td determined by kTf, <t< kT^ + Td with k = 0, 1, ... during which the coordinates possess a definite 
sign, that is they are either non-negative or non-positive depending on the initial conditions. (We recall that Aq is 
non-negative by definition.) At the beginning of each interval the values are w.l.o.g. given by Aq(kTi,) — Ago = 
and Aq{kTii) = Ago > 0. Positive (negative) velocities (7„(fcTf,) > {q^iikTb) < 0) with qmax{kTb) > qmin{kTb) > 
{(JminikTb) < qmax{kTb) < 0) result in non- negative (non-positive) amplitudes gn(i) > (qn{t) < 0) during intervals 
kTb < t < kTb + Td. That is, all oscillators are a.t t — kTb situated at the position q„ — 0, corresponding to the 
minimum position of the on-site potential, and Aq(kTb) and Aq{kTb) attains its maximum and minimum respectively. 
During kTb < ^ < kTb + Td/2, the quantity Aq{t) monotonically decreases resulting at i = kTb + Td/2 in Aq{t) = 
while the monotonically increasing quantity Aq{t) reaches its maximum. During kTb + Tb/2 < t < kTb + Td both Aq(t) 
and Aq{t) < monotonically decrease attaining at the end of the interval kTb + Td their minima Aq{kTb + Td) — 
and AqlkTb + Td) = -Aq{kTb). 

For the subsequent stage of length Tb — Td, when the amplitudes qn{t) have opposite sign compared to the previous 
interval, the motion of Aq{t) and Aq(t) starts with the same values as at the beginning of the previous interval, viz. 
Ago = and Ago > and the oscillator at the lattice site that previously supported qmax (and (jmax) possesses now 
the minimal amplitude qmin (and minimal velocity qmin) and vice versa. However, Aq(t) and Ag(i) resemble the 
behaviour of their counterparts during the previous interval. 

In order to establish conditions for the non-existence of breather solutions we consider the behaviour of Aq{t) and 
Ag(i) with initial values Ag(0) — and Ag(0) > w.l.o.g. on intervals 

h ■.kTb<t< kTb + Td, with fc = 0, 1, 2, .... (9) 

during which the amplitudes qn{t) are either exclusively non- negative or non-positive (see above). Ag(i) is smooth 
on the intervals Ik- Furthermore, on each interval Ik it holds that 

Aqit^ (kTb + Td/2)) = Aqi-t- (kTb + Td/2)) (10) 

Aqit- {kTb + Td/2)) = -Aq{-t- [kTb + Td/2)) (11) 

with fc = 0, 1, ..., viz. Aqit) is even with respect to tk = kTb + Td/2 whereas Ag(i) is odd. 

Further, exploiting the symmetry features and periodicity of Aq{t) and Ag(i) one obtains the following relations: 

Ag((fc + l)Tb) = Aq{kTb + Td) = Aq{kTb), (12) 

Aq{(k + l)Tb) = -Aq{kTb + Td) = Aq{kTb) (13) 

with Ag(0) = and Ag(0) > 0. Crucially, the relations ()12p and ()13p constitute necessary conditions to be satisfied by 
breather solutions. Thus, for given values Aq{kTb), Aq{kTb) at the beginning of intervals Ik the solution Aq{kTb + Td), 
Aq{kTb -\- Td) at the end of intervals Ik can be utilised to derive a first recurrence (Poincare) map for which breathers 
constitute fixed points. 

The time evolution of the difference variable Aq{t) is determined by the following equation 



<fAq 



= - [U'iqmax) - U'{q„iin)] 

'~f\Qriiax ~^ Qmin) ~r J^\f^)\Qmax "~ Qmin) 



/ ^ \y^max \_^ \Q7nax-\-j Qmax) f V \Qniax Qniax—j)\ 

^rain Y*^ \Qmin-\-j Qinin ) i V yq^nin Qmin—j)\j ■ V-'-^j 

Notice that the direct driving field A{t) has no impact on Ag(i). Regarding the maintenance of localisation the 
inequality Aq{t) > constitutes a necessary condition. Regarding the equal sign, for localised solutions, such as 
breathers, where the oscillators perform in-phase motion (and/or out-of-phase motion) Aq{t) is zero only at instants 



of time when the oscillators, whilst performing periodic motion inside their potential wells, attain simultaneously 
the position q^ — ai the minimum of the on-site potential. Conversely, if A.q{t) approaches zero in finite time no 
localised pattern persists at all. 

For the forthcoming derivations of estimates we facilitate the following statement: 

Lemma: For soft potentials U{q) with two inflection points qi^± consider the interval 

Is ;= [qi,qr] with gi,_ < qi, and, qr < (?»,+ . (15) 

Then it holds that for any pair x,y £ Is with x < y 

[U'{y) - U'{x)] >5s{y~x)>0 (16) 

where the constant (5^ > is given by 

Ss = min U"{q) = min [U"{qi), C/"(g.)] . (17) 

For hard potentials consider the interval Ih '■— [xi^Xr] with — oo < xi, Xr < cw. Then it holds that for any pair 
x,y G Ih with x < y 

[U'{y)-U'{x)]>6h{y-x)>0 (18) 

and the constant Sh > is given by 

Sh = mmU"{q) = U"{0). (19) 

q&lh 



Proof: Consider the expression 



y~x 

By assumptions ([5]) and ([7]) we have that on intervals Is and Ih it holds that U'{y) > U'{x) for y > x. Therefore the 
expression F[x, y) is positive. Furthermore, by virtue of the mean value theorem there exist a point z in (x, y) such 



that 

U'{y) - U'ix) 



U"{z) > min U"{q). (21) 



y — x qeisJh 

One has for soft potentials mingg/^(t/"(g)) = min [U"{qi), U"{qr)] > 0, and therefore it holds that 

U'{y) - U'ix) > min [U"{qi), U"{qr)] {y - x) ^ 6siy - x) > 0. (22) 

Similarly for hard potentials by the assumption ([2|) one has inmq^j^{U"{q)) — U"{0) > 0, so that 

U'{y) - U'ix) > C/"(0)(y - x) ^ 5h{y - x) > 0. 

completing the proof. 

D 

Remark: To apply Lemma to the case of soft potentials with a single inflection point qi > one proceeds along the 
lines given above for Lemma 1 considering the interval {—oo, qr] and qr < qi- The positive constant 6s is given by 
,5, = min,e/^ [[/"((?)] = t/"(<z,). 

In the following we present conditions for which Ag(t), associated with a breather solution satisfying the conditions 
listed above, exponentially decays which rules out the existence of breather solutions to Eq. ([T]). 

Theorem: Let the relation (7/2)^ > Wq — Bmax > be valid with cjq = 6s and ujq = 6h for soft and hard on-site 
potentials given in Eq. (IT71) and P^ respectively. Then it holds that Eq. ^ does not possesses breather solutions. 

Proof: We prove the assertion by contradiction. That is we suppose that Eq. ([T|) exhibits breather solutions associated 
with periodic functions Aq{t + T^) = Aq{t) and and Aq{t) — Aq{t + Ti,) satisfying the necessary conditions in p^ 



and (|13p . Using the conditions in ([2]) and (HI together with the Lemma enables us to bound the r.h.s. of Eq. (|14l) on 
each of the intervals /j,, fc = 0, 1, .., defined in (jH]), from above as follows: 






\U [Qniax ) ^ \Qriiin )\ "T \ Qrnax Q'rnin ) 
-^V J \Qrnax Qrnin) 
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^max 


j=i 


\ 



V'{qraax+j " qynax) + V (q. 



max Hmax—j) 



>0 



/ 



* \Qmin-\- j Qriiin) i V [Qmin Qniin—j } 



V 

dAq 



<o 



<o 



max ^min) ^[Qniax Qjnin } 
2 



7nax Hm,i7i J 



-(W(, - B,nax)Aq 



7- 



dt 



Therefore Ag(t) is bounded from above by the solution of 

d'^Aq 



df' 



(Wg - Brnax)Aq - 7 



dAq 
dt 



(23) 



The solution to Eq. (|23p is given for (7/2)^ > ojq — B^ax > on each interval J^. by 



Aqk{t) 



exp 



(-!') 



Ago, A; + ^Ago,; 



W^ 



■ sinh(W^t) 



Ago.fe cosh(M^i)] , 



(24) 



and 



Aqk{t) = exp(-^t) 



2VF 



Ago^fc 



sinh(T4^t) + A4o,fc cosh(H/i) 



(25) 



where the index k refers to the interval I^ and W = y^{j/2)'^ — {ujq — Bmax)- Due to the Eqs. ([12]) . (|T3)) . fulfilled by 
breather solutions, the following recursion relations are true 



Ago,fe+i = Agrj. (fcTft + Td) 
Aqo.k+i = -Aqk (kTb + Td) 



(26) 
(27) 



with starting values A(7o.fc=o = and Ago.fc=o > 0. Using the latter recursions and the notation Qj — Aq{jTi,) and 
Pj = Aq(jTii) with j — 0, 1, ... we cast the solution in form of a first recurrence (Poincare) map 



Pj+i 



M 



P. 



where the matrix M is given by 



with entries 



M = E 






(c + ^s) 



E 



exp(-2^<i 



C = cosh(W^Trf) 
Sj = sinh(l^Td). 



(28) 

(29) 
(30) 
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FIG. 1: Time evolution of Aq{t) exhibiting on average exponential decay in accordance with inequaUty (|24[) for the system 
given in Eq. (|32p . The values of the parameters are 7 = 0.35, Aq = Bo = 0.99, flA ~ Sis — 2.0, Qa ~ Q% = 0, and k — 0.5. 



For the determinant of M one obtains 



det M = -E^ {C^ - 5^) = -E^ 



(31) 



As I det M\ < 1 the Poincarc map is contractive and as a consequence Aq{t) and Aq{t) exponentially decay approaching 
zero which is in contradiction to the condition of periodic behaviour of non-vanishing Aq{t) = Aq(t + T^) and 
Aq{t) — Aq{t + Ti,) and the proof is complete. 

D 

Conclusively, as Aq converges to zero our theorem provides conditions that rule out the existence and/or formation 
of breather solutions. This statement holds not only for single breathers but also for multi-breathers, which can 
be regarded as arrays of single breathers, where the oscillators perform in-phase or out-of-phase periodic motion 
with respect to a reference oscillator [3,l2g| provided Aq{t) involves the same two oscillator. As proved in [23] the 
only available multibreather solution in one-dimensional lattices are those with relative phase (in-phase) and tt 
(out-of-phase) between the lattice sites. 

Remarkably, the process of exponential decay takes place regardless of the amplitude of the external field A{t). 
Furthermore, exponential decay happens for any kind of attractive interaction potential V{u). As far as hard on-site 
potential U{q) is concerned, only its curvature at the bottom, U"{0), plays a role for the decay process and the larger 
is the curvature the faster is the exponential decay while increasing the amplitude of the parametric driving B^ax 
has the opposite effect. Note that in order that the theorem applies the latter has to fulfill the constraint Bmax < i^q- 
Importantly, the result holds for general driving fields. Interestingly, in our upper bound the decay rate turns out to 
be independent of the initial distribution of the amplitudes {g„(0)}. They influence the amplitude of the decay of 
Aq{t) though. 

We stress that the hypothesis (7/2)^ > luq — Bmax > can be satisfied for arbitrarily small values of the damping 
strength 7 as for given Wg = t/"(0) for hard on-site potentials (wq = ^^^[U"{qi),U"{qr)] for soft on-site potentials) 
the amplitude of the parametric driving field Bmax can be tuned to control the infimum of 7 complying with the 
inequality. Hence, 7 can be sufficiently small compared to a characteristic frequency of the system (which is e.g. 
given by oscillations near the bottom of a potential well with frequency determined by U"[Q)) so that the system's 
dynamics is kept away from the overdamped limit. 

For an illustration of the exponential decay of an initially localised solution we choose for the hard on-site potential 

The interaction potential is harmonic and is given by 

V{qn - qn-l) = -^{qn - qn-lf, 

where the interaction radius is taken as r = 1 amounting to linear nearest-neighbour interaction and the coupling 
strength is uniform, i.e. k„ = k. As the external fields are concerned we consider time-periodically varying fields and 
set for the direct driving field A(t) 



A{t) = AQS\n{Q.At + Q\) 




FIG. 2: Spatio-temporal pattern of the coordinates qn{t) for a chain consisting of A'' 
couphng. The parameter values are as given in Fig. [T] 



100 oscillators with nearest-neighbour 



with amplitude ^Oj frequency Ha and phase Q'^. Similarly, for the parametric driving field B{t) we choose 

B{t)^Bosm{^Bt + e%), 

with amplitude Bq, frequency fls and phase Q'g. 
The corresponding lattice system is given by 

Qn = -qn ~ ql + K (qn+i ~ 2qn + q„-i) - ^Qn 
+ Aosm{nAt + QA)+BQam{nBt + e%)qn. 



(32) 



In our simulation the system comprises N = 100 oscillators and periodic boundary conditions are imposed. We plot 
in Fig. [T]the temporal behaviour of Aq{t) for the system (|5^ starting from a localised single hump solution peaked 
around the site n = 50 associated with initial conditions g,i(0) — l/cosh(n — 50)) and qn{0) — 0. Ag(i) exponentially 
decays on average which is in accordance with the inequality (j24p bounding the amplitude of A.q{t) from above. The 
corresponding spatio-temporal evolution of qn{t) is shown in Fig. [2] corroborating the exponential decay of a spatial 
pattern. Eventually the dynamics settles on a spatially uniform state and the oscillators perform identical motion, 
qn{t) — q{t) and quit) = q(t) for all n, entailing that the oscillators decouple. The oscillators undergo periodic motion 
on a limit cycle supporting periodic oscillations of the variables qn{t). 

In conclusion, we have studied the persistence and formation of non-homogeneous patterns represented by breather 
solutions in general nonlinear damped and driven one-dimensional lattice systems. Sufficient conditions, in terms of 
the values of the parameters, have been provided which assure that no non-uniform state can exist. To be precise, it 
has been proven that the difference between the maximal and minimal ampUtudes of the lattice oscillators of a non- 
uniform state decays exponentially fast. In this way we have proven that creation and/or preservation of spatially 
(localised) patterns is impossible. Notably our results are independent of the number of oscillators. Conversely, 
rigorous quantitative conditions are identified under which discrete breathers can exist in general driven and damped 
lattices at all. Furthermore, our generic results on the non-existence of localisation and pattern formation in general 
nonlinear lattice systems open the way for further research on the emergent dynamical scenarios, in particular features 
of synchronisation, in coupled cell networks. Given that we have provided quantitative criteria in parameter space for 
the existence/nonexistence of discrete breathers the current work is also expected to stimulate further experimental 
studies of breathers in nonlinear damped and driven lattice systems. 
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